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SUMMARY 

A semi empirical analysis of the equation for incompressible fluc- 
tuations in a turbulent fluid, using similarity relations for round 
subsonic jets with uniform exit velocity, is used to predict the shape 
of the time-averaged fluctuation-pressure distribution along the mean- 
velocity boundary of jets. The predicted distribution is independent 
of distance downstream of the nozzle exit along the mixing region, in- 
versely proportional to the distance downstream along the region of 
mean-velocity self-preservation, and proportional to the inverse square 
of the distance downstream along the fully developed region. 

Experimental results were in fair agreement with the theory. How- 
ever, the measured fluctuation-pressure distributions were found to be 
very sensitive to changes in jet temperature and jet-nozzle profile, 
especially near the nozzle. These factors are not included in the the- 
ory. Increased jet temperatures produce increased pressure fluctuations 
and violation of similarity conditions. Nozzle-profile modifications 
may lead to violation of the uniform-exit-velocity requirement imposed 
in the theory. 


INTRODUCTION 

Some previous theoretical studies of noise from round subsonic jets 
(refs. 1 to 3) have related the radiated acoustic power to the strength 
distribution of noise sources along a jet. Measurements to confirm the 
theories have not been performed. On the other hand, measurements of 
pressure fluctuations along the mean- velocity boundary of a jet have 
been performed (ref. 4) ; but no attempt has been made to predict theo- 
retically the shape of the fluctuation-pressure distribution. This lat- 
ter problem appears to be as important as the former in evaluating the 
noise characteristics of jets and is considered herein. The previous 
theoretical studies in references 3 and 5 are particularly relevant to 
the following discussion. 
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The author is greatly indebted to Professor Herbert S- Ribner for 
suggesting the possibility of a mathematical attack on the present prob- 
lem and for his many enlightening discussions of jet noise theory. 


THEORY 


Consider a continuous, perfect, homogeneous, inviscid medium contain- 
ing a region of turbulence possessing a mean motion, but free of sources 
of matter and external forces. The appropriate exact equation for fluc- 
tuations which occur within the medium is, in rectangular Cartesian co- 
ordinates, x-^(i = 1, 2, 3), 


c^p 

Sx? 


bt 2 = " 6x i^ x j 


(p v i v j) 


( 1 ) 


where p is pressure, p is density, vp are velocity components , and. 
t is time. (All symbols are defined in the appendix.) 


A pressure-sensitive microphone placed along the mean-velocity bound- 
ary of a jet will respond to temporal fluctuations of the pressure asso- 
ciated with incompressibility, as well as to temporal fluctuations of 
the pressure associated with compressibility. The former fluctuations, 
called "pseudosound," are not radiated, whereas the latter fluctuations 
(sound) are propagated acoustically. Because the pressure fluctuations 
very near the jet are dominantly associated vith incompressibility 
(refs. 1 and 6), the detected fluctuations siould be described essenti- 
ally by the incompressible part of the inhomc geneous wave equation (l), 
namely 
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where the density p Q is constant over space' and time. Set p = p + p 
and v k = v k + v k , where p and v k are temporal-mean values and p 
and v, are fluctuations associated with incompressibility. (The flue- 
tuations p and v^ correspond, respectively, to pseudosound and tur- 
bulence.) Prom the resulting equation, subtract its temporal mean. The 
result is 
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(2b) 
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Equations (2) possess the form of Poisson's equation. Hence, the solu- 
tion of equation (2b) is 


p(x; t) = - 



d 3 y 



(4) 


where x denotes the observation point, as shown in figure 1, y denotes 
source points, the volume integration is over all space (but for prac- 
tical purposes, only over the source region), and the velocity components 
are functions of y and t. By performing a double partial integration, 
equation (4) becomes (ref. 7) 


p(x;t) 


PO 

4rt 



*ij(y;t) 


(jj - x i )(y 1 - Xj) 



d 3 y 


(5) 


Equation (5) expresses the instantaneous pressure fluctuation at one 
observation point (xp,X 2 ,X 3 ) as a function of simultaneous instantaneous 
velocity distributions throughout the turbulent region. However, in 
order to associate p with measurable velocity distributions, the veloc- 
ities, hence p, must be expressed in the form of time averages. The 
^^ e ~ avera ^ e P i- s zer °j but, in general, the space-time covariance 
of p is not zero if the spatial separation, or time delay, or both, 
are not excessively large. Thus, consider 

p(xjt)p(x’ ; t 1 ) 



(y i -x i )(y i -X .) 


y- x 


d3 y 
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d 3 y’ 


The vectors x,x’ denote different observation points; y,y’ denote dif- 
ferent source points; and t,t' denote different times. In the present 
instance interest is confined to the mean- square-pressure fluctuation. 
Hence, x' = x, and t' = t. Consequently, 

p 2 (x; t) 
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where j(y 1 ;t) represents a space covariance. Usually the 

space covariance of a quantity is expressed as a function of location in 
space and the separation of the observation points (source points in the 
present instance). Thus, consider the transformation 


y = y 


y f = y + z 


where z connects source points (see fig. l) . Then, equation (6) 
becomes 


p 2 (3J;t) 


2 

p 0 

16*' 
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'l' i j(y;t)Ti» kI (y+ zjt) 


(yi- - x ,i) (yk + z k - x kXn + z z - x z) 




|y+ z - 5?| : 


(7) 


The integrations indicated in equation (7) appear considerably more 
complex than those performed in reference 5 because any product 

y' ;t) is a function of y and z, rather than z alone. 
Performing the integration does not appear presently feasible. This 
prevents the computation of absolute values of p 2 . However, a semi- 
empirical analysis of equation (7) can be performed in a manner analo- 
gous to that performed by Lighthill (ref. 8) for the case of acoustic 
power radiated from a region of turbulence. This should permit a deter- 
mination of the shape of the fluctuation-pressure distribution along the 
jet. In Lighthill' s case the entire turbulent region was characterized 
by a single characteristic velocity, frequency, and length. Moreover, 
the entire turbulent region contributed to pressure fluctuations at any 
given far-field observation point. In the present instance the situation 
is somewhat different, but analogous. Whereas in Lighthill 's analysis a 
single value of a single quantity (total acoustic power) was sought, in 
the present instance a distribution of values of a single quantity (time- 
averaged pressure fluctuations) is sought. Thus, the entire turbulent 
region cannot be characterized by a single set of values of character- 
istic quantities. Rather, the characteristics of each region of the flow 
must be specified separately. Moreover, although the entire region of 
turbulence contributes to p at any near-field observation point, the 
observed pressure fluctuations will be very strongly determined by tur- 
bulence in the flow immediately adjacent (see eq. J/7) and ref. 7). This 
results because of the very rapid attenuation of p 2 as a function of * 
the source-point - observation-point separation |y - x|, or |y + z - x| . 
Possibly the only related experimental data are presented in figure 16 
of reference 4. These data show that for any given frequency the tur- 
bulence and pressure fluctuations peak at adjacent points along the jet 
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length. Thus , if it is assumed that the pressure fluctuations p at a 
given observation point along the jet boundary are determined essentially 
by the turbulence within an adjacent limited volume x of the jet, the 
analogy with Lighthill's analysis is complete. 

Assume that the characteristic volume t is proportional to an aver- 
age eddy volume at the downstream location of the observation point. The 
eddy volume is approx imately equal to the volume over which 

^i j (y (y + / 0, where y denotes a fixed point and 1? denotes 

a movable point. Let R represent a characteristic separation of the 
fixed^point and the observation point; R corresponds to |y - x| and 
|y + _z - x| in equation (7). (For purposes of the analysis the effect 
of z on the values of other quantities may be disregarded, since only 
characteristic values are of interest.) Because the region influencing 
p(x;t) extends only over the volume x, it will be assumed that the in- 
tegrations over ^ and z in equation (7) need only extend over x, or 
possibly over a volume proportional to x. The remaining factors 

yi - x i , yt + z i - x i 

~T~ 777 and TZ Z — represent directionality factors which are 

|y - x| |y + t - x| 

to be disregarded in an analysis of the present type. 

The foregoing argument leads to the following expression for p2 
along the jet boundary: 

p 2 * p oVu t2r ' 6 

It is usual to express the velocity components in a round jet in cy- 
lindrical (rather than rectangular Cartesian) coordinates r, cp, z. Let 
Zl> t)z* ^3 correspond to r, cp, z, respectively (see fig. l). Let u^, 

u 2> u 3 represent, respectively, the radial, azimuthal (dimensionless), 
and axial components of velocity in the cylindrical coordinate system. 
Velocity components in the two coordinate systems are typically related 
by 


dyi 

^ Uk 


where u, tj in the cylindrical system correspond to v, y, respectively, 
in the rectangular system. The dimensional forms of the velocity com- 
ponents in the two coordinate systems differ only by trigonometric fac- 
tors which are to be disregarded in the present analysis. 


When is expanded in terms of the velocity components Vj_ 

(see eq. (3)), and hence u^ (see preceding paragraph), second-, third-, 
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and fourth-order covariances of the turbulent-velocity components re- 
sult. The magnitudes of these covariances are generally not well known 
(ref. 5 and ref. 9, ch. VIII). However, because the mean- velocity com- 
ponent u 3 is generally very much larger than any of the other co mpo- 
nents, the term i 3i = ujUf, and hence the covariance m u 3 > 

should be expected to dominate. Thus, the dimensional representation of 

p^ may be rewritten as 


~2 2 _ 2 ~ - ■ ■ 2_-0 

P a Po u 3 u i u Z TlR 


( 8 ) 


Proportionality (8) applies at any given point along the jet bound- 
ary. The variation of g2 as a function of £ 3 is sought along the 
jet boundary, where | 3 is the jet axial coordinate of the observation 
point measured from the nozzle-exit plane. The situation is analogous 
to that of applying Lighthill's acoustic power relation (ref. 8) for 
similar jets. In addition to the usual requirements that the flows be 
geometrically similar and that the mean-flow and turbulent-flow Reynolds 
numbers be large, the present analysis implies that the eddies distrib- 
uted along the jet be geometrically similar. 

Along the jet, t <x and R a (ref. 10) . Hence, 


c~ ~ 

P - 


( 9 ) 


For round subsonic jets where the mean- v cloc ity distribution over 
the nozzle exit is uniform, similarity relations have been found (ref. 
11, pp. 183-186) which relate velocity components up, u± to the geom- 
etry of the jet. The relations and their associated regions of validity 
are as follows: 


Region 1 : 1 t} 3 / d <5 4 

^3 = t ) t ) 

Region II : 8 < d < 20 to 50 



(10a) 


(10b) 


/ uw 
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Region III: 20 to 50 <S rij/d. 



- u o % ® a P \Ty 

^3 x 7 


(10c) 


where Uq is exit velocity; d is nozzle-exit diameter; 

where (ru) is the virtual origin of the jet; and a,p 
^ 0 


1 * ^3 " ^ 3^0 
1,3. Corre- 


sponding relations for the other components have not been determined. 
However, up « U3, and ug = 0. The magnitudes of the unlisted turbu- 
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lence covariances (containing the azimuthal components) are likely not 
greater than those given in equations (10). 


In the notation of the preceding paragraph, region I is the mixing 
region. Region II is a region of self-preservation of the mean-velocity 
distribution. Region III is a region of self-preservation of the mean- 
and turbulent- velocity distributions called the region of fully devel- 
oped turbulence. 

The distribution of p 2 as a function of r) 3 is obtained by sub- 
stituting the similarity relations (10) in proportionality (9). Thus, 

Region 1 : 1 < rij/d < 4 

p 2 /( 

Region II: 8 < < 20 to 50 

p 2 / (po”o) 

Region III : 20 to 50 ^ n 3 /d 

p 2 / (po^o) 

where () denotes some form of space average over the volume x. 

Although F, G, f, and g are functions of tj 3 , the dependence is 
very weak compared with the dependence on t^. The primary dependence 
of p 2 on T13 is contained in the coefficients of the velocity func- 
tions. Also, by virtue of the assumption that the observed pressure 
fluctuations at a given distance downstream of the nozzle exit are 


<&o) ' ( f 2<5 3^ (lla) 

I 

- (i ) 2 < f2 °^> (llb) 

■ (tT < f V) (llc) 
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determined by turbulence at the same distance downstream of the nozzle 
exit, the geometrical coordinate of the source point may be re- 

placed by the coordinate £3 denoting the observation point. It follows 
that, to a first approximation, 

Region 1 : 1 < £ 3 / d < 4 


F/^o 4 ) * (^/ d )° ( 12a ) 

Region II: 8 < £^/d < 20 to 50 

¥/(&$ * (^/ d )" 2 ( i2b ) 

Region III: 20 to 50 < |^/d 

P 2 /(pQ U o) a ^3/ d ) _i ( 12c ) 

Thus, the theoretical distribution of the time-averaged pressure 
fluctuations along the boundary of the mixing region (l <; 13 /d <; 4) is 

flat, similar to that of the acoustic- source strength (ref. l), that is, 
the power radiated per unit length of jet. However, along the boundary 
of the fully developed turbulent region (20 to 50 £ l^/d), the pressure 

fluctuations decay considerably less rapidly than does the acoustic- 
source strength , the decay law being proportional to £ 3 ^ for mean- 
square pressure and proportional to ^ for source strength. 

EXPERIMENT 

Distributions of time-averaged fluctuating pressures along the mean- 
velocity boundary of round, cold, subsonic air jets were measured using 
3- and 5- inch-diameter nozzles. Profiles of the nozzles are shown in 
figure 2. The measured pres sure -fluctuation distributions are shown in 
figure 3. The distribution along a turbojet -engine exhaust (ref. 4) is 
also shown. The engine nozzle consisted of a straight tailpipe termin- 
ated by a spherical convergent iris. The engine exhaust pressure- 
fluctuation distribution is for a hot jet (T ** 1070° F) at a slightly 
supercritical value of nozzle total-pressure ratio (p/pq 2.2). A 

sound-level recorder possessing a quasi-peak averaging characteristic 
was used for the engine tests, whereas a meter having a full-wave 
rectified averaging characteristic was used for the air- jet tests. 
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The distributions in figure 3 have been replotted in the appropriate 
dimensionless form in figure 4. Corresponding theoretical curves are 
also shown. The distribution for the engine exhaust has been reduced by 
3 decibels to account for the difference between the meter averaging 
characteristics. The shapes of the theoretical and experimental distri- 
butions tend to agree. However, the levels of the engine- exhaust and 
air-jet distributions differ, and the distribution for the 5-inch- 
diameter nozzle does not obey the (^3/d) 0 law along the mixing region. 

The level difference may result partly from inaccurate correction 
for the difference between the meter averaging characteristics and partly 
from dissimilarity of the microphone locations in the different tests. 

(Note that p 2 is a strong function of near the nozzle (ref. 4).) 

However, dynamical and geometrical dissimilarity of the jets caused by 
the jet temperature difference and the nozzle-profile difference may have 
a greater influence on the level difference than the preceding factors. 
With regard to the effect of jet temperature, it should be noted that 
entropy fluctuations were neglected in the mathematical analysis and 
hence are neglected in the dimensionless ordinate used in figure 4. 

Also, jet temperature influences the space distribution of p. The ef- 
fect of jet temperature on the jet geometry is shown in reference 12. 

An increase in jet temperature tends to shorten the length of the jet 
core, increase the rate of jet spread, and increase turbulence levels. 

(in this connection, the microphone traverse angle with respect to the 
jet axis was the same (10°) in all tests.) 

Some results are shown in figure 5 which illustrate the effect of 
different values of jet temperature on the fluctuation-pressure distri- 
bution. These tests were for the same nozzle (5-in- diam. , fig. 2) , so 
that nozzle geometry was not influential. The microphone traverse angle 
was also the same (11.5°) for all tests. Results are shown for a sub- 
critical value of nozzle pressure ratio (fig. 5(a), P/pq = 1.31) and 
for a slightly supercritical value (fig. 5(b), p/pq - 1.99). As shown, 
for example by equations (12), the appropriate dimensionless ordinate 

for the subcritical-pressure-ratio case is p 2 /pqU^. In the present 

instance (fig- 5 (a)) this ordinate representation has been applied in 
the modified form 


* SPL 2 - SPL-l + 20 log 
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where SPL is sound pressure level (referred to 2X10 -4 dyne/cm^), the 
subscript 1 is associated with the lesser value of jet temperature, 
and the subscript 2 is associated with the greater value of jet temper- 
ature. By virtue of the representation the effect of the 

jet-velocity difference on pressure-fluctuation distributions measured 
for constant nozzle-pressure ratio should be eliminated. The representa- 
tion j/P/p^ is "based on the assumption of incompressible flow. For 

supercritical values of pressure ratio the flow is certainly not incom- 
pressible. For nonisentropic jets operating; at supercritical values of 
nozzle-pressure ratio, the pressure ratio j/pq has been found to be a 

correlating parameter for acoustic power gereration at constant jet tem- 
perature (ref. 13). Applying this parameter in the present instance (fig. 
5(b)), the appropriate ordinate representation for supercritical values 

of nozzle-pressure ratio is |/p^/P, which in figure 5(b) has been ap- 
plied in a modified form analogous to the subcritical ordinate 
representation. 

The fluctuation-pressure levels plotted in figure 5 were repeatable 
within 1 decibel, except possibly for > 50 inches, where atmospheric 
winds were likely to influence the jet flow considerably. Although the 
temperature difference was small (82° F), the fluctuation pressures for 
the warmer jet were consistently higher near the nozzle than those for 
the cooler jet. Downstream of the nozzle, where jet temperatures tended 
to equalize, the pressure fluctuations (dimensionless form) tended to eq- 
ualize. Thus, a definite effect of jet temperature on near-field pres- 
sure fluctuations is indicated. 


The failure of the fluctuation-pressure distribution for the 5- inch- 
diameter nozzle to obey the (| 3 /d)° decay relation illustrates the ef- 
fect of the nozzle profile on the fluctuation-pressure distribution. It 
should be recalled that the velocity- similarity relations (10) apply only 
if the mean-velocity distribution over the rozzle exit is uniform. This 
condition is certainly not satisfied in the case of the 5-inch-diameter 
nozzle. Schlieren photographs of the flow from the 5- inch-diameter noz- 
zle have shown the emergent boundary layer to be highly turbulent. In 
contrast, schlieren photographs of the flow from the 3- inch-diameter 
nozzle have shown the emergent boundary layer to be quasi-laminar down- 
stream to a vena contracta, from which plane the jet expands an d the 
mixing region becomes highly turbulent. In this case the mean- velocity 
distribution over the nozzle exit can be expected to be fairly uniform. 
The convergence of the nozzle produces a negative pressure gradient which 
tends to preserve laminar flow in the boundary layer. Sketches of the 
expected flows from the two nozzles are shown in figure 6. On the dimen- 


sionless basis 


(VP/pqUo); 


pressure fluctuations of larger magnitude 
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should be expected near the exit of the 5-inch-diameter nozzle. Also, 
the ratio of the maximum averaged fluctuation pressure to the value of 
the averaged fluctuation pressure at the nozzle-exit plane should be ex- 
pected to be greater for the 3-inch-diameter nozzle by virtue of the fact 
that the emergent boundary layer has transformed from quasi laminar to 
turbulent. These expectations are confirmed in figure 4. The effect of 
dissimilarity of the nozzle profiles on near-field pressure fluctuations 
is considerable. 


CONCLUSIONS 

The incompressible-fluid theory for near-field pressure fluctuations 
in conjunction with jet similarity relations appears to account for the 
general shape of the distribution of time-averaged pressure fluctuations 
along the mean- velocity boundary of a round subsonic jet with uniform 
exit velocity. However, jet temperature, which is not considered in the 
theory, significantly affects the values of the fluctuations and the 
fluctuation distribution, especially near the nozzle. Higher jet temper- 
atures result in larger pressure fluctuations. The jet-nozzle profile 
can also significantly affect the fluctuation-pressure distribution, es- 
pecially near the nozzle, by introducing a nonuniform exit-velocity 
distribution. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 22, 1960 
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APPENDIX - SYMBOLS 
jet-nozzle-exit diameter 
volume elements 

velocity distribution functions 

self -preserving velocity distribution functions 

stagnation pressure 

static pressure 

ambient pressure 

characteristic separation of source point and observa- 
tion point 

temperature 

time 

jet-exit (core) velocity 

velocity components in cylindrical coordinate system 

velocity components in rectangular Cartesian coordinate 
system 

rectangular Cartesian coordinates of observation point 
(origin at center of nozzle exit; see fig. l) 

rectangular Cartesian coordinates of source point 
(origin at center of nozzle exit; see fig. l) 

rectangular Cartesian coorlinates of source point 
(origin at another source point) 

cylindrical coordinates of source point (origin at 
center of nozzle exit; s?e fig- 1) 

cylindrical coordinates of observation point (origin 
at center of nozzle exit; see fig. l) 

dens ity 


t? 

i 

c 

c 


ambient density 
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T 

+ij 

0 


characteristic local eddy volume 

zVj + - Vj 

space average over volume x 


Subscripts : 

a^p = 1,3 components in cylindrical coordinate system (l^ radial 

component; 3, jet axial component) 

i, j,k,2 * 1,2,3 components in coordinate systems (in rectangular 

Cartesian coordinates: 1,2 , components normal to 

jet axis ; 3, jet axial component. In cylindrical 
coordinates: 1, radial component } 2, azimuthal com- 

ponent, 3, jet axial component.) 


Superscripts : 

time average 

incompressible fluctuation 

' source point or observation point generally different 

than that associated with corresponding unprimed 
quantity 

vector 
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(b) Nozzle-pressure ratio, P/Pq, 1*99 (supercritical). 

Figure 5. - Effect of Jet temperature on fluctuation- pres sure distribution along mean-velocity boundary, 
Nozzle-exit diameter, 5 inches; ambient temperature, 71° F. 








